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Abstract 

We compute the threshold uncertainties due to unknown masses of the Higgs 
bosons on the predictions for the intermediate and unification scales, Mj and Mjj 
respectively in SO(10) models. We focus on models with separate breaking scales for 
Parity and SU(2)r symmetries since they provide a natural realization of the see- 
saw mechanism for neutrino masses. For the two step symmetry breaking chains 
, where left-right symmetric gauge groups appear at the intermediate scale, we find 
that parity invariance of the theory at the unification scale drastically reduces the 
GUT threshold effects in some cases. Including the effects of the intermediate scale 
thresholds ,we compute the uncertainty in the above mass scales and study their 
implications for proton lifetime and neutrino masses. An important outcome of our 
analysis is that if the currently favored nonadiabatic MSW solution to the solar 
neutrino puzzle is accepted , it will rule out the SU{2) L X SU(2) r XU(1)b-lX SU(3) c 
as an intermediate symmetry for SO(10) breaking whereas the intermediate symmetry 
SU(2) L XSU(2) R XSU(A) C , is quite consistent with it. 
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I. Introduction: 

The grandunified theories 1 (GUTs)provide an elegant extension of physics be- 
yond the standard model. The requirement that the gauge couplings constants in these 
theories become equal at the GUT scale (My) lends them a predictive power which 
makes it possible to test them in experiments such as those looking for the decay of 
the proton. The most predictive such theory is the minimal SU(5) model of Georgi and 
Glashow 1 , where the SU(5) symmetry breaks in one step to the standard model. The 
only new mass scale in this model is the My which can be determined by the unifi- 
cation requirement using the low energy values of any two gauge couplings from the 
standard model. One then predicts not only My, but also the remaining low energy 
gauge coupling constant (say sin 2 8w)-H is well known that for the minimal SU(5) 
model, they lead to predictions for the proton life-time as well as sin 2 9w both of which 
are inconsistent with experiments. 

This however does not invalidate the idea of grandunification and attention has 
rightly been focussed on SO(10) 2 GUT models which can accomodate more than 
one new mass scale. Supersymmetric SU(5) 3 models also belong to this class. In 
this class of two mass scale thoeories,the values of low energy gauge coupling con- 
stants can determine both the mass scales again making these theories experimentally 
testable. The determination of the values of the new mass scales become more precise 
as the low energy values of the gauge coupling constants become better known. It is 
therefore not surprising that the recent high precision measurement of a str ong and 
sin 2 9w at LEP 4 once again revived interest in grandunified theories 5 . 

Supersymmetric SU(5) theories have been studied with the goal of predicting the 
scale of supersymmetry breaking 5 . These models however do not have any room for a 
nonzero neutrino mass nor natural generation of adequate baryon asymmetry, whereas, 
the SO(10) model is the minimal GUT scheme that provides a frame-work for a proper 
understanding both these problems. In this paper,we concentrate on the SO(10) mod- 
els with a two step breaking to the standard model and study the threshold effects 
on the predictions for the two new mass scales i.e. My and M/.In order to appreciate 
the significance of our work, it is worth pointing out that in SO(10) models the scale 
My as usual is related to proton decay whereas the intermediate scale is related to 



2 



neutrino masses if the intermediate symmetry is either of the left-right symmetric 
6 groups SU(2) L XSU(2) R XG C where G c is SU{A) C or SU(3) c XU(1) B -l-^ the neu- 
trino mass is determined independently (say, from the MSW solution to the solar 
neutrino puzzle), then the see-saw mechanism determines the range of the required 
intermediate mass scale. The viability of a given SO(10) model will then depend on 
both the value of Mj obtained from renormalization group analysis as well as the 
uncertainties in this value arising from threshold corrections. 

Let us discuss the kind of SO (10) models we will study here. As is wellknown,the 
SO(10) group contains the maximal subgroup SU(2) L XSU(2) R XSU(A) c xD, D being 
a Z 2 symmetry which implements the parity transformation^ well as particle-anti- 
particle transformation). We will refer to this symmetry as D-parity. The actual 
nature of SO(10) model depends on what symmetry appears at the intermediate 
scales (i.e. between the GUT scale and M w ).The most interesting SO(10) models 
are the ones where the symmetry breaking to the standard model occurs in two 
steps,with either of the left-right symmetric groups SU(2) L XSU(2) R XSU(A) C or 
SU(2) l XSU{2) r XU(1) B -lXSU(3) c (denoted henceforth by G 224 andG 2213 respec- 
tively) as the only intermediate symmetry. These are also the theories for which def- 
inite predictions can be made. Our work will focus on them. Note the absence of 
D-parity at the intermediate scales .Use of Higgs multiplets belonging to 45 and 210 
representations to break SO(10) can lead to such a scenerio as was pointed out in a se- 
ries of papers in 1984 by us in collaboration with D.Chang 7 . Let us briefly recapitulate 
some other motivation for considering such models. 

One of the attractive features of the SO(10) models is that they provide a natural 
understanding of the neutrino masses via the see-saw mechanism 8 . It has however 
been noted that, the see-saw mass matrix does not follow naturally in models where 
D-parity and SU(2) R breaking scales (M R ) are same. On the other hand, if the D- 
parity breaking scale M P is such that M P 3> M R then the see-saw formula emerges 
naturally 9 . This is perhaps the most compelling motivation for requiring the D- 
parity breaking scale to be significantly larger than Mr. There are however other 
motivations from cosmology.If M P =M R , there arise domain walls bounded by strings 
at the epoch when SU{2) R XD symmetry breaks down .They dominate the mass 
density of the universe making it hard to understand the successes of the big bang 
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picture. Such problems do not arise in SO(10) models with separate D-parity breaking 
scenerios. Furthermore, exact D-parity leads to n B = rig. In SO(10) models where the 
baryon asymmetry of the universe arises from Higgs boson decays, the ratio ns/n 7 
receives an additional suppression (Mp/Mu) 2 on top of its small value predicted in 
generic GUT mo dels. This mechanism would prefer scenerios with D-parity breaking 
scale to be above that of SU (2) R and at the GUT scale . We are of course fully aware 
that, if baryon asymmetry arises from the decay of heavy Majorana neutrinos,the 
above constraint does not apply. 

A complete two-loop analysis of the predictions for sin 2 9w and proton life-time 
in this class of SO(10) models was carried out in ref. 10. Depending on the nature of the 
Higgs boson spectrum used to implement the symmetry breaking and the nature of 
the intermediate symmetry groups, the intermediate scales and the associated physical 
implications were discussed. The mean values of the mass scales will be taken from 
this paper. 

A basic limitation of all grandunified theories is that all mass scale predictions 
are subject to uncertainties arising from Higgs boson thresholds 11 . It has therefore 
been argued 12 that since the Higgs bosons in question belong to large representations 
in SO(10) theories, one might worry that the mass scale predictions derived from two- 
loop calculations are completely unreliable. In other words, even if sin 2 9 w and a str0 ng 
are very precisely known,the unification scale M v and the intermediate scale Mj will 
have large uncertainties. It was however subsequently pointed out that this need not 
always be true;for instance, if an SO(10) model has an intermediate symmetry group 
SU(2) L XSU(2) R XSU(A) C XD, (G 224D ) the GUT threshold uncertainties in sin 2 9 w 
exactly cancel out 13 . This result of course holds only if the intermediate symmetry 
is G224D and does not apply to the more interesting models with separate D-parity 
breaking ;it also does not say anything about the uncertainties due to intermediate 
scale thresholds.lt does however give rise to the hope that existence of symmetries 
may reduce the net impact of threshold uncertainties. In any case if the grandunified 
theories are to be useful, threshold effects must be calculated. In this paper ,we begin 
this program for the two SO(10) theories and hope to extend it to other models later 
on. The main results of this paper have already been reported earlier 14 . 
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In ref.14 and in the present paper,we adopt the following approach. Using 
the evolution equation for the coupling constants, we express the My and Mj in 
terms of the known low energy parameters a str ong,sin 2 9w and a em and the threshold 
corrections due to Higgs bosons .Since the LEP results have considerably reduced 
the experimental uncertainty in sin 2 6w as well as a s ,the main uncertainty comes 
from the arbitrariness associated with the Higgs boson masses and the theoretical 
uncertainties in the scales My and Mj can be computed. The final magnitude of 
the uncertainty depends o n how far the scalr masses are split from the symmetry 
breaking scale. Using the standard model as a guide,we assume that constraints of 
one-loop radiative corrections and unitarity bound on tree level amplitudes would 
allow the scalar boson masses to be a factor of 10 on either side of the symmetry 
breaking scale. We also present results for a wider splitting of (30) ±x for illustration 
even though we believe this to be rather unlikely. 

We consider the symmetry breaking chains: 

(A) SO(iOHG 22 4^G std ; 

(B) SO(10HG 221 3^G std . 

We compute the Higgs boson threshold effects on the uncertainty in intermediate 

+ 2.7 

scale [Mc/Mq] in case (A) to be 10- 1 - 4 and that in the value of grandunification scale 
i.e. in {My /My) to be 10- 1 7 .This corresponds to a maximum value for t p in case(a) to 
be 10 40 years for a s =.ll and 10 38,4 years for a s —.l. For case (B),we find the uncertainty 
in (M R /M° R ) to be 10-^ and that in (My/My) to be at most 10 ±2 .Note that the 
threshold uncertainties are much less than the estimates of ref.6.For a s =.ll,we obtain 
an upper limit on t p in this case to be 5 x 10 36 years. 

Using our results in combination with the see-saw formula for neutrino masses, we 
find that the presently favored nonadiabatic MSW solution to the solar neutrino puz- 
zle rules out the SO(10) model (B) which has SU(2) L XSU(2) R XU(1) B - L XSU(3) C as 
an intermediate symmetry.In our opinion this is an important result since this will be 
the second GUT model that is being definitively ruled out by experiment. The sym- 
metry breaking chain (model (A)) is however quite consistent with data. This paper is 
organized as follows: in sec.2,we present a derivation of the formulae for the threshold 
uncertainties for the model (A);in sec.3,we derive the contribution of the various Higgs 

5 



multiplets to these uncertainties and give our estimate of these effects. Sec. 4, uses the 
results of sec. 2 and 3 to derive the same results for the model(B);in sec. 5, we derive 
the implications of our results for the solar neutrino puzzle and in sec. 6, we discuss 
the implications of our results for proton decay and the effect of adding extra Higgs 
multiplets on our result. In sec. 7, we close with some concluding remarks. 

II. The Formula for the Threshold Corrections : model (A): 

Let us now proceed to derive the equations for the threshold uncertainties. We 
will illustrate the technique for the model(A). We start with the standard renormal- 
ization group equations for the evolution of the gauge coupling constants: 

/ daA aiaj 1 2 

In eq.(l),a!i = |^ and the one loop coefficient a; = — + |ra s + where n g is 
the number of fermion generations and T(s) is the contribution of the Higgs bosons. 
The bij are the two loop coefficients, which are not needed here. At each symmetry 
breaking threshold,we use the following maching conditions 14 : (we assume that the 
group Gi breaks to the group Gi at the scale Mr) 

Oi(Mj) oij 12vr 

In eq.(2), \{ = Tr(9Y) 2 +Tr(9f ) 2 Zn(M i/ /M / ); QY are the generators of the lower 
symmetry Gi for the representations in which the heavy Gauge bosons appear;^ is 
the same for the superheavy Higgs bosons. 

Let us now apply the formulae in eq.(l) and (2) to the SO(10) models described 
earlier as (A) .Note that since the D-parity has been broken in both cases at the GUT 
scale, the Higgs multiplets needed to implement the symmetry breaking in model (A) 
are in 210,126 andlO dimensional representations of the SO(10) group. We denote by 
Mjj,Mc and Mz the three symmetry breaking scales and they arise from the vev's of 
the above three Higgs multiplets respectively. The D-parity breaking manifests itself 
in the mass of the submultiplet of the 126 representation A L (transforming as (3,1,10) 
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under G224) being different from the submultiplet A R transforming as (1, 3, 10) under 
the same group. We postpone any further discussion of the Higgs bosons to the next 
section. Let us now derive the formula for the threshold corrections. 

Using eq(l) and (2) and the standard Georgi-Quinn- Weinberg type analysis,we 
find the following expressions for sin 2 6 w and a s : 



16- (V - = AJn (^) + Auln {^f) + F c + F x + T s (3) 



where 

A c = (8a 3 — 3a 2 L — 5ay — 604 + Sa 2R + 3a 2L ); 

A v = (6a 4 - 3a' 2R - 3a 2L ); 

F c = 3lnf3(a' 2L - a 2L + a 2R - a' m - 2a 4c + a 4c ); 

Fx = _ (V 3 )(^3c - §A 2L - § A y + §Aj£ - |A^ - fA^) 



™*f . , fl 3\_ DI M c , D M v 



stn 2 6 w - - = £ c Zn-/ + + G c + G A + T e (4) 

a em V 8/ M z M z 

where 

£ c = (5(a 2 L - ay) - (5a 2L - 2a 4c - 3a 2i? )); 
B v = (5a 2L -2a 4c -3a 2i? ); 

G c = ln(3(5a' 2L — 3a 2R — 2a 4c — 5a 2L + 3a 2ij , + 2a 4c ); 

r» _ 5/\C \C ,3\U , 2\[/ \[/ \ 
— gl^y - ^2L + 5 A 2i? + 5 A 4c _ A 2L) 

In the above expressions, a^a^ and a" denote the the evolution coefficient for 
the gauge couplings between M w toM c ,M c to M P and M P to M v respectively. M P 
is the scale of D-parity breaking. The T's denote the two loop contributions, which 
do not contribute to threshold uncertainties to the leading order and will therefore 
be omitted henceforth. The values of the a^'s for model (A) are: a' 2L = a' 2R = a' 2R = 
11/3, a 2l = —3, a' 4c = —23/3, a 4c = —14/3, a Y = 41/10, a 3c = —7; 

They will be used in the numerical estimates of the various effects. Using eq(3) 
and (4), we can express the mass scales Mq and My in terms of the low energy 
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parameters and the threshold contributions. The uncertainties in the low energy pa- 
rameters are experimental and can be estimated to be small as we show. The threshold 
contributions buried in the A's intrduce the theoretical uncertainty having to do with 
the fact that the heavy Higgs masses are unknown. 

Let us first address the uncertainties due to the experimental errors in sin 2 9w 
and a s which we take 16 as: 

sin 2 6 w = .2334 ± .0008; a s = 0.115 ± .007 (5) 
Denoting by C = ^(^f - |) and C, = ^(sin 2 w - §), we get 



im M °\ a / CqBu — C\Au 

(Aln— W = A( ) (6) 



(\i A/ CoB c - CiA c 

(Aln— W = M BuAc _ AvB ) (7) 

Using eq.(3) and(4), and the values of various a^'s given above to evaluate the 
A's and B's , we can now estimate the uncertainty in the quantities and |pf where 
the quantities with subscript denotes the values corresponding to the mean values of 
sin 2 9w and a s : 



2rt = 10± - 025 ; ^ = 10± - 22 (8) 

We will find that these uncertainties are small compared to those arising from 
unknown masses of heavy higgs bosons, thanks to precision experiments from LEP 
e + e~ collider. 

Again using eq.(3) and (4) and the expressions for A,B,F and G, we can derive 
the following expressions for the uncertainties in M\j and Mc arising from threshold 
effects only: 



AlU O = f P M+f U M+f C M-(A M /A e ) (f p d + f u 9 + f c e )+(AlnM c ) expt (a s , sin 2 6 w ) 



Aln (H) = ~T e ifeP + ff>U + feC) + ( AlnM uU P Ms, stn 2 6 w ) (9) 



where 



and 



f M = ( i - T'r ) inp (loa) 



^ M {s(a' 2L ^\' 2R -2a' 4c )) {m 



; M l,18(a' 2i + a' 2ii -2agj ^ 



f 6 - ~ { (a> L + a> R -2a> c ) ) ^ 



j 9 ~{ 6(a' 2L + a' 2R -2a> 4c ) ) ^ ] 



f c = I p2L« - a' 2R ) + (5/3)A^(a 2L - a' 4c ) + Xg(a' 2R + (2/3)a 4c - (5/3)a 2L ) \ 
9 6 \ a' 2L + a' 2R - 2a' 4c J " 

(10/) 

/5/3a y + a 2L - 8/3a 3c \ 
= 1 - 7-— —, (10^) 

A e = 5/8 ( i — r- ) £ (Wh) 

\a 2L + a 2R ^a 4C y 

where 

B = ((3/5a 2iJ + 2/5a4 C - a 2L )(5/3a y + a 2L - 8/3a 3c ) - (a 2L + a' 2R - 2a' 4C ){a Y - a 2L )) 

In eq(10) ,/3 = Mu/Mq- if the intermediate symmetry is G224D,one can still use 
the above general expressions, after dropping the ff and terms in eq(8) and(9).(In 
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this case of course, all a = a" .) Note that, in this case f$ and f[ vanish. To see that 
this is what one expects from the results of Parida and Patra(ref.l3), we note that 
in their work (ref.l3),the uncertainty in Mq was assumed to be zero. Using this and 
bringing back the sin 2 9w and a s terms to the equation (9), fg and ff terms can be 
identified as the GUT threshold uncertainty in sin 2 9w and therefore their vanishing 
in the G 2 2ad limit was what was established in ref.12. 

III. Survival Hypothesis and Estmation of the Threshold Uncertain- 
ties: 

In order to give a numerical estimate of these uncertainties,we need to know the 
masses of the physical Higgs bosons;more specifically, what submultiplets are at what 
mass scale. This can be done using the survival hypothesis for the Higgs bosons 17 . The 
basic assumption of the survival hypothesis is that only a minimal number of fine 
tunings of the parameters in the Higgs potential are done as required to ensure the 
hierarchy of the various gauge boson masses. In the case at hand we need to fine tune 
only two parameters since we have only a two step breaking. The survival hypothesis 
then says 17 that a submultiplet of the Higgs multiplet of the GUT group,that acquires 
a vev to break a given subgroup Gi is stuck at the symmetry breaking scale. The other 
submultiplets which transform as complete irreducible representations under Gi get 
pushed to the next higher scale. Using this,we find the scales for the Higgs boson 
masses . 

For the case of model(A),the Higgs multiplets needed are 210,126 and two 10- 
dimensional multiplets. They are responsible for the three symmetry breaking scales 
M[/, Mj = Mq and M\y respectively.As is well-known, the 210 multiplet also breaks 
the D-parity symmetry. Using the survival hypothesis, the scales of the different Higgs 
multiplets can be obtained and they are listed in table I and II. In table I, we list the 
Higgs bosons with masses around My; in Table II, the Higgs bosons with masses near 
Mj are listed. 

The U-submultiplet is the Goldstone mode corresponding to the superheavy 
Gauge bosons and are omitted in computing the threshold uncertainties. 

Let us now give their contributions to the various A^'s and to final uncertain- 
ties.Defining rji = (Inj^), we find: 

10 



A 2 l = \ U 2 R = 6 + 30% + 3077s + 20t7 C (1 1) 

Here,we have assumed that M Sl = M Sfl by left-right symmetry. Similarly, d and 
(2 effects are combined and denoted by C-The coefficients in front of the different 77's 
are simply the Dynkin indices of the different multiplets under the different gauge 
groups. For instance,for t]2l, the Dynkin index is that of the gauge group SU{2) etc. 



\ U 4C = 4 + 2 VH + 2 Vs + 32r] Co + 24% + 24r? c + 4% (12) 

Repeating the similar procedure for the intermediate scale,we find the contri- 
bution to A c 's to be as follows: 



Ay- = 1/5 (3r/ + 2 VRl + 4% 2 + 16% 3 + 32% 4 + 64% 5 + 24% 6 + 14) (13) 



= V* (14) 

^3c = 1 + VRi + 5% 2 + 5% 3 + % 4 + 5% s (15) 
The threshold contributions to Aln(jf-) and Aln(jf-) can now be written down 



AZn(— -) = (1.29Zra/3 + .0259 (-2 + 2r] w + 4% 6 + 2r/ 210 ) + .0017 (65 - 15^ + 769%)) 
Mz 

(16) 



Aln(— -) = -min(3 - .O62I7710 - .124?7 126 - .062r7 210 + .03% - .496% (17) 

Mz 

In deriving eq(16) and(17),we have assumed that the Higgs multiplets belong- 
ing to a single SO(10) supermultiplet have the same mass at a given scale 18 . This 
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implies that r]^ = rjs = t\yi% , = Vc an d a U Vr, are equal. While there can be 
deviations from this degeneracy assumption, their contributions will not change our 
results noticeably. Secondly, we allow the e^'s to be between (l/10)and (10) as well as 
1/30 and 30. The first case corresponds to allowing the Higgs self-scalar couplings 
to range from 10~ 2 to 10 2 times the gauge coupling and is better motivated by the 
analogy of standard model than the second choice, although, we present the results for 
both cases. As mentioned above the value of M& L is always kept lower than M Vl which 
corresponds to taking (5 less than one. Using eq(16) and(17), we compute the uncer- 
tainties in the intermediate mass scales as well as the Mf/.We present these results in 
table III. There are two possibilities;one when the uncertaintyin Mc is maximized and 
another when the uncertainty in Mjj is maximized. We find the maximal uncertainty 
in jjtj- to be 10- 1 7 whereas that in |j£ to be 10- 1S from Higgs boson threshold effects 
alone. We will study the impact of our results on the predictions of the SO(10) model 
in sec. 5. 

IV.Threshold Corrections for Model (B): 

Let us now turn our attention to discussing mo del (B), where the SO (10) symme- 
try first breaks down at the scale M v to SU{2) L XSU{2) R XU{l) B - L XSU{3) c which 
subsequently breaks down at scale M R to the standard model. Again as before the D- 
parity symmetry is broken at the GUT scale. The Higgs multiplets necessary to imple- 
ment this chain are 45,54,126 and 10 dimensional ones. The D- Bodd component of the 
45-dim. Higgs breaks the GUT symmetry.lt has been pointed out that 19 , without the 
presence of 54,45 will break SO(10) down to SU(5)XU(1) rather than G2231 .However 
as far as the threshold corrections are concerned,in the limit of exact degeneracy 
the contribution of the submultiplets of 54 exactly cancel. Let us now present the 
equations for the threshold corrections to Mjj and M R in this case. Using the same 
notation as in eq.(8) and (9), (except that we replace M c by M R in eq(9)),we give the 
expressions for the various f 's and A's below. Defining A'jj = 8a 3c — 3a' 2L — 3a' 2R — 2a' BL 
and B' v = 5a' 2L — 3a' 2R — 2a' BL ,we get 

1 5 
A 6 = -n-(8a 3c - 3a 2 L - 5a Y - -fy-(a 2 L - a Y )) (18a) 
A u B u 
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A M = \ \ - -^(8a 3c - 3a 2L - 5a Y )j (186) 



ff = Infi \J^( Sa 3c - 6a 2L - 2a BL) - J^( a 2L ~ Q Bl)J ( 18c ) 



f¥ = ^( 8X l - 6A& ~ 2Ag L ) + ^r{\ U 2L - \ U BL ) (18d) 



= e^(8A£ - 3A* - 5Ap) + A-(A« - Ay) ; (18e) 



fu = W [l ~ ^H&4 - 6a' 2 ' L - 2a' Bi j (18/) 



/m = g^(8A£ - 6A& - 2Ag L ) (18«?) 



/£ = -£(8A£-3A&-5A£) (18/i) 

In order to evaluate the threshold corrections, we need the values of gauge cou- 
pling evolution coefficients a^'s as well the mass scales of the physical Higgs bosons. 
We have a 2L = a m = a' 2R = -7/3,A"bl = 7 a' 3c = a 3 ' c = -7,a' 2L = -3,a BL = 11/2 
a 3c = -7;a 2L = -19/6, a Y = 41/10. 

Next we will use the survival hypothesis to determine the mass scales of the var- 
ious Higgs bosons. In table IV,we give the Higgs bosons with masses of order ofM^.We 
do not include the components of 54-dim multiplet since in the degenerate multiplet 
approximation , their effects cancel out exactly. The Higgs fields with masses of order 
M R are only two in number and are therefore listed in the text. 



There are only two Higge multiplets at scale M R ; they are A^ + (l, 2^3/5, 1) and 



0(2, — |y3/5, 1) where the numbers within the bracket refer to their transformation 
property under the standard model. Their contribution to the various A's are given 
below: 
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A^(126) = 2A VHlL + 12r] H2L + 6 VH3 + 677^ + 16^ 5 + 2^ 6 (19a) 

A&(126) = A^ L (126) 

A£(126) = 30r] Hl + 677^ + ^]h 3 + 4^ 4 + 24^ 5 + 77^ + r] Hs (196) 

A^ L (126) = 12 VHl + 677^ + I677J/3 + 16t7h 4 + ?7ff 7 + ^ 8 (19c) 

A^(10) = A^(10) = (19d) 

AL(10)=A^(10) = ^+i»i) (19e) 

A&(45)=2i7 a (19/) 

A^(45) = 27753 (19») 

A 3 l/ C (45) = 3^ (19/i) 

A^(45) = (19j) 

A? = 6/5 + 24/577^ + ^ (19fc) 



A* = ^ (19/) 

\* = o (19m) 
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Using eq(18) and (19) and the values of aj's,we get for the threshold contribution 
to the uncertainties in My and Mr the following expressions: 



Aln(^) = .0685/n/3 - .I71r] 126 - .049^ - .039r] w - .177^ - .146?^ (20) 
Mz 



Mr 

AZn(— = .095/n/3 - .083r/ 126 + .033^5 + -062?? 10 - .06^ + .22r] Rl (21) 
Mz 

In order to evaluate the possible uncertainties, as before we keep /3 less than one 
and allow e v to vary between .1 to 10 in one case and 1/30 to 30 in the second case. The 
results are given in table III. WE see that in this case the threshold uncertainties are 
much less than in model (A). The maximal uncertainty in ^ is 10-- 3 whereas that 



in is 10 ±2 .This has 



R 



We also wish to note at this point the uncertainties in My and Mr arising from 
the errors in a s and sin 2 9 w : Using the same formula as in eq (6) and (7), we find for 
model (B): 

^% = lO^ 25 ; ^ = lO^ 18 (22) 
M v ' M R y ! 

V.Solar Neutrino Puzzle and SO(10): 

In this section, we study the implications of the results derived in this paper 
for solar neutrino puzzle. As is well known,one of the most interesting resolutions of 
the solar neutrino deficit is the so called MSW matter oscillation mechanism 20 . In 
this mechanism,resonant enhancement of the oscillation of v e to either or v T takes 
place in the solar core for a range of values of the Am 2 and the mixing angle #.In the 
so-called high mass (adiabatic) solution, the value of Am 2 is of order 10~ 4 et> 2 with 
sin 2 9 ~ 0.02 — .6 whereas in the nonadiabatic solution, we instead have Am 2 sin 2 29 ~ 
4X10~ 8 eU 2 with Am 2 ~ 10~ 6 eU 2 - 8 x 10^ 8 eU 2 . The combination of Chlorine 21 and 
Kamiokande 22 and initial Gallium data 23 seems to point towards the nonadiabatic 
solution 20 . Either of the cases seem to fit quite well with the see-saw picture for the 
neutrino masses in the SO(10) model 24 provided one assumes D-parity breaking 9 . In 
the presence of D-parity breaking,the hierarchical quadratic mass formula for neutrino 
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masses follows naturally ( a fact, which appears not to have been well appreciated by 
many theorists). The point,briefly is that, due to the presence of couplings of Higgs 
bosons of type 126. 126. 10. 10, a vev of A R induces a vev of the 126 submultiplet A L of 
order 1 ^ M - .This leads to a direct mass for all left-handed neutrinos of the same order 
invalidating the conventional see-saw mechanism formulae. If however ,the D-parity is 

2 

broken at the GUT scale,the A L vev becomes only of order 9 m ^X R which is smaller 

M u 

than the see-saw contribution to the neutrino masses. 

After the radiative corrections are taken into account 25 , the formulae for neutrino 
masses are: (assuming generation mixings to be small) 

77? 2 



rnu, = (-07) ^ (236) 

M N 

m VT = (.18)^ (23c) 
M N 

In order to find the neutrino masses,we need to know which is given by 
Mn = (f/g)M R .lt can be argued on the basis of vacuum stability 26 that / < g.The 
mean value of M R has been obtained from two loop analysis of the two SO(10) models 
in ref.10 and 27. For case (A) we have Mq ~ 10 11 - 5 GeV,whereas for model(B),we have 
M R ~ 10 9 GeV.The uncertainty in the exponent due to the error in a s and sin 2 6w 
is about ±.025 in model(A) and ±.18 in model(B). Including this and the threshold 
uncertainties, we find the minimum value of neutrino masses (corresponding to f=g ) 
to range between the following values: 

Model(A):m^ = 6 x 10- 12 eV - 2 x K)- 7 eV ;m Vlt = 4 x 10- 7 eV - K)- 2 eV; 
m Vr = .006eV - 180eV. 

Model(B):m^ = lO' 7 eV - 10~ 6 eV = lQ~ 2 eV - Iti^eV and m Vr = lOOeV- 
lkeV. 

We see that if the adiabatic solution is ruled out as is currently believed, then 
model(B) will be ruled out by the solar neutrino experiments and only model(A) will 
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be acceptable. This is an important result in our opinion since there are no other 
uncertainties one can hide behind to save this model. 

VI. Higgs Boson Related Uncertainty in Proton Decay and Stability 
of the Threshold calculations: 

In this section ,we discuss two questions: i) the uncertainty in the predictions for 
proton decay and ii) the effect of adding extra Higgs bosons to a GUT theory on 
the above calculations. First, we discuss the predictions for proton in the two SO(10) 
models under discussion. Again using the results of ref.10 and 27 we find that for 
a s = .11, the value of My = 10 15 ' 8 GeV. The uncertainty in a s , leads to an uncertainty 
of order 10 ± 22 (see eq(8) and (22)) multiplying the above value. We predict the proton 
life-time for the model(A) to be r p = 1.6 x 10 35± - 7± - 9 -6-s years.For the model(B),we 
find, t p = 1.6 x io 35± - 7±L± - 8 years. 

Let us now turn to the question of the stability of our results. It is sometimes 
stated that any additional Higgs multiplet added to a GUT model will add to the 
already existing uncertainty. However, in a recent paper 28 , it has been shown by one 
of the authors that if the additional Higgs multiplet does have a vev or has vev in a 
gauge direction which has been broken by a Higgs field with the same representation 
content, then threshold effects from such multiplets always cancel in sin 2 9w or the 
intermediate scales. This lends a degree of stability to the above calculations. 

VII. Conclusions: 

In conclusion,we have presented a detailed analysis of the threshold effects due 
to the unknown masses of the Higgs bosons and shown their effect on the numer- 
ical predictions for the values of the unification and the intermediate scales in two 
SO (10) models with a two step breaking. To the best of our knowledge, this is the 
first time such an analysis has been carried out for the SO (10) models. An interesting 
outcome of this analysis is that the nonadiabatic MSW solution to the solar neu- 
trino puzzle is inconsistent with the model (B) which has an intermediate symmetry 
SU(2) L XSU(2) R XU(1) B - L XSU(3) C . 
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Table Caption 



Table I. The Higgs Bosons at Mass scale Mjj. 

Table II. The Higgs bosons with masses at Mc\ The numbers within the paren- 
thesis refer to the representation content under SU(2) L XU(1) Y XSU(3) C The multi- 
plet ^arises from the 0(2, 2, 0) and the R-multiplets arise from the multiplet A#(l, 3, 10) 

Table III. In this table, we present our results for the threshold uncertainties in 
the intermediate scale and the unification scale for different values of e^.The first four 
lines correspond to the case where the uncertainty in Mj is maximized whereas the 
last four lines correspond to the case where the uncertainty in 

Table IV. Higgs bosons with masses of order Mu . 

Table I 



SO (10) Representation 


G 2 24 sub multiplet 


10 


H(l,l,6) 


126 


Co(2,2,15),S(l,l,6),A L (3,l,10) 


210 


S L (3,1,15),S R (1,3,15), 
Ci(2,2,10),C 2 (2,2,fO) 
C 3 (1,1,15),S'(1,1,1,),U(2,2,6) 
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Table II 



SO(10) representation 


G213 submultiplet at Mq 


10 


0(2,-^/3/51,1) 


126 


^(1,^/3/51,3), 
i? 2 (l,^3/5|,6 
^(1,-^/3/51,6) 
^4(1,-^3/51,3) 
^5(1,^3/51,6) 
i? 6 (l, 2^3/5,1) 



Table III 



Symmetry Breaking chain 


M H /Mu or Mh 1 Mi 


M 7 /M 7 ° 




C224 


1/30 to 30 


+4 

10-21 


+1.2 

10-2.5 


C2213 




+ .9 
10-4 


+ .1 
10-2 


C224 


1/10 to 10 


+2.7 
10-14 


+ .8 
10-17 


C2213 




10^3 




C224 


1/30 to30 


+4.2 
10-2.2 


+ 1.2 
10-2.5 


C2213 




+ .5 
10-2 


io ± - 2 


C224 


1/10 to 10 


+2.8 
10-1.8 


+ .8 

10-17 


C2213 




+ .3 

lO-o 


io ± - 2 
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Table IV 



SO(10) representation 


G2213 content of the heavy boson 


10 


71(1,1,1^/3/2,3) 
T 2 (l, 1,-^3/2, 3) 


126 


H lL (3, 1,-173/2, 6) 
tf 1R (l,3,+!73/2,6) 

/f 2L (3, 1,-|7 3 A 3) 
i? 2jR (l,3,+| v / 3/2,3) 

iJ3(2,2,-f73/2,3) 

F 4 (2,2,+f73/2,3) 
# 5 (2,2,0,8),/i 6 (2,2,0,l) 
H 7 (l, l,|73/2, 3) 
if 8 (l, l,-|7 3 /2, 3) 


45 


^(1,1, 0,8), ^(3,1,0,1), ^(1,3,0,1) 
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